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Abstract 


The  induced  noise  control  parameter  is  defined  in  terms  of  the 
ratio  of  the  stored  energy  in  a  master  dynamic  system,  when  it 
is  coupled  to  an  adjunct  dynamic  system,  to  that  stored  energy 
when  the  coupling  is  absent.  An  induced  noise  control  parameter 
comprises  of  two  factors.  The  first  is  defined  in  terms  of  the 
ratio  of  the  indigenous  loss  factor  of  the  master  dynamic  system 
to  the  virtual  loss  factor.  The  virtual  loss  factor  is  the  sum 
of  the  indigenous  loss  factor  of  the  master  dynamic  system  and 
the  induced  loss  factor.  The  induced  loss  factor  measures  the 
influence  of  an  adjunct  dynamic  system  that  is  coupled  to  the 
master  dynamic  system.  The  first  factor  is,  then,  by  definition 
less  than  unity.  The  second  factor  is  defined  in  terms  of  the 
ratio  of  the  external  input  power  into  the  master  oscillator, 
when  coupled,  to  the  external  input  power  in  the  absence  of  this 
coupling.  This  second  factor,  which  is  positive  definite,  may 
either  exceed  unity,  be  unity,  or  be  less  than  unity.  Both 
factors  in  the  induced  noise  control  parameter,  however,  are 
critically  dependent  on  the  global  coupling  strength.  The 
global  coupling  strength  is  again,  a  ratio.  This  ratio  is  that 
of  the  energy  stored  in  the  adjunct  dynamic  system  and  in  the 
coupling  to  the  energy  stored  in  the  master  dynamic  system.  The 
investigation  of  this  ratio  and  its  relationship  to  the  induced 
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loss  factor  is  explained  and  graphically  illustrated.  The  modal 
coupling  strength  is  determined  via  the  global  coupling 
strength.  In  the  statistical  energy  analysis  (SEA)  the  value  of 
the  modal  coupling  strength,  by  definition,  lies  below  unity. 

In  the  energy  analysis  (EA)  herein  developed,  the  modal  coupling 
strength  may  exceed  unity.  This  excess  occurs  in  the  energy 
analysis  (EA)  when  the  coupling  is  strong  and  the  damping 
assigned  to  the  adjunct  dynamic  system  is  low.  The  damping  in 
the  adjunct  dynamic  system  is  conveniently  defined  in  terms  of 
an  associated  modal  overlap  parameter.  To  reconcile  the 
statistical  energy  analysis  (SEA)  with  the  energy  analysis  (EA) 
in  focus,  the  associated  modal  overlap  parameter  in  the  adjunct 
dynamic  system  must  necessarily  exceed  a  threshold. 
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I.  Introduction 

Damping  treatment  is  often  employed  in  order  to  diminish  the 
response  of  an  externally  force-driven  reverberant  dynamic 
system.  The  statement  is  linked  to  the  definition  of  damping  in 
terms  of  the  loss  factor  //ol®)'  namely 

?7<,W=[n“W/ffl£°W]  ;  <(<»)=[n°(0)/<»)7„W]  ,  (d 

where  Ylgi^)  is  the  power  input  from  the  external  force-drive, 

E°{o))  is  the  stored  energy  in  the  dynamic  system  and  {co)  is  the 

frequency  variable  [1-3] .  The  stored  energy  in  a  dynamic  system 
is  an  eligible  measure  of  its  (quadratic  average)  response;  see 
Fig.  1.  From  Eq.  (1)  it  follows  that  increasing  the  damping, 
which  entails  an  increase  in  the  loss  factor,  may  result  in  a 
decrease  in  the  response. 

Utilizing  Eq.  (1),  a  stored  energy  reduction  scheme  may  be 
proposed:  The  dynamic  system  in  focus  -  the  master  dynamic 
system  -  is  modified;  e.g.,  by  appropriately  coupling  it  to 
another;  see  Fig.  2.  The  adjunct  dynamic  system,  to  which  the 
master  dynamic  system  is  to  be  coupled,  is  suitably  designed  to 
change  the  loss  factor  as  perceived  by  the  external  force-drive. 
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from  rio(a})  to  rj^  (fi>)  with  [rio(‘>)iV.  (g?)]<1.  The  change,  by  definition 
and  in  the  vein  of  Eq.  (1),  yields 

/7vW=[ne(0)/ftj£o(6>)];  Eo{(o)=\^e{(o)l (Ori^{(o)[  ,  (2) 

where  the  corresponding  changes  in  the  stored  energy  and  in  the 
external  power  are  from  Eq{(o)  to  from  YV^ip^  to  > 

respectively  [1-10] .  The  external  force-drive  in  this 

consideration  is  assumed  to  remain  unchanged  under  the 
modification  to  the  dynamic  system  as  just  proposed.  The  induced 
noise  control  parameter;  designated  ^(<»),  may  then  be  expressed, 
from  Eqs.  (1)  and  (2),  in  the  form 

[£„(<»)/ £;(<»)] ;  #„W=[/7„W/?7vW]  -(3) 

where  the  external  input  power  ratio  7r{o})  is  defined 

^(®)=[n.(®)/n°(®)]  ■ 

Employing  Eqs.  (3)  and  (4)  one  may  speculate  that  if  the  external 
input  power  is  uninfluenced  by  the  modification;  i.e.,  if 
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7r{co)=i 


,  (5a. 1) 


and  if,  in  addition,  the  modification  successfully  achieves  a 
loss  factor  ?]y{co)  such  that 

the  induced  noise  control  is  highly  beneficial;  namely,  the  noise 
control  parameter  ^{o))  is  small  compared  with  unity;  ^(®)«1.  The 

condition  stated  in  Eq.  (5a. 1),  however,  may  not  hold.  After 
all,  a  change  in  the  loss  factor,  perceived  by  an  external  force- 
drive  in  a  dynamic  system,  may,  under  certain  conditions,  be 
expected  to  change  the  external  input  power.  Then,  increasing 
the  degree  of  damping,  in  the  manner  here  prescribed  and  stated 
in  Eq.  (5b),  may  result  in  an  input  power  ratio  7r{(o)  which  may 
either  stay  as  stated  in  Eq.  (5a. 1),  exceed  unity  or  fall  below 
unity 

7r{co)>l  ,  (5a. 2) 

or 

7r{(o)<l  .  (5a.  3) 

In  the  absence  of  heavy  handedness  in  the  construction  of  an 
adjunct  dynamic  system  7l{(a)  usually  exceeds  unity;  i.e.,  7T{(o)  is 
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usually  coiranensurate  with  Eq.  (5a. 2).  Then  the  highly  beneficial 
noise  control  that  Eqs.  (5a. 1)  and  (5b)  are  promising,  is 
mollified  as 


where  use  is  made  of  Eq.  (3)  [11].  The  Principle  that  underlines 

this  mollifying  influence  -  the  mitigator  of  many  a  noise  control 
panacea  —  is  attributed  to  Chatelier  [12] .  Le  Chatelier 
Principle  in  Noise  Control  Engineering  is  largely  the  subject 
matter  of  this  report. 

Focusing  attention  on  Fig.  2,  the  relationships  that  are 

governed  by  the  conservation  of  energy  are  employed  to  define  a 
number  of  loss  factors,  and  the  global  and  modal  coupling 

strengths.  With  the  appropriate  interpretation  of  these 

definitions  one  may  derive  an  estimation  of  the  external  input 
power  ratio  /r(fy)  and,  hence,  provide  a  more  definitive  estimate 
for  this  quantity  than  Eq.  (5c)  does  [11-13].  Then  the  noise 

control  parameter  stated  in  Eq.  (3),  may  be  properly 

estimated.  ,  Largely  these  estimates  are  cast  in  terms  of  mean- 

values,  as  prescribed  by  Skudrzyk  [14] .  Computational  illustra¬ 
tions  are  cited. 
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II.  Relationships  stemming  from  the  conservation  of  energy 

Focusing  attention  on  Fig.  2,  one  may  supplement  Eq.  (2)  with 
the  following  definitions  derived  off  the  conservation  of  energy 
(power) .  This  conservation  may  be  expressed  in  the  forms 

rio  (<»)  =  Vo  Eo  («)] ;  n,  {co)  =  Vi  («)[®  (")]  =  Vs  {oi)[(o  E,  («)] ; 

ne(<y)=no(6?)+n,(<y);  ne(6>)=7v(®)[6>^oH]  /(^a) 

and  hence 

TiM=Vo(<»)+ri,{m),  w  =  [i + to  W/  (41’ 

3o(®)  =  [-Ei®)/-E„(®)]  -(6b) 

where  is  the  portion  of  the  external  input  power  dissipated 

in  the  master  dynamic  system,  and  n^(o)f  E^ico)  and 

portion  of  the  external  input  power  dissipated,  the  energy  stored 
and  the  loss  factor  in  the  adjunct  dynamic  system,  respectively, 

and,  finally,  is  the  ratio  of  the  energy  stored  in  the 

adjunct  dynamic  system  to  that  in  the  master  dynamic  system  [1- 
3] .  (The  coupling  elements  are  considered  as  part  of  the  adjunct 
dynamic  system;  e.g.,  the  stored  energy  E^ip)  includes  the  stored 
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energy  that  resides  in  the  coupling  elements.)  The  stored  energy 
ratio  3o(<y)  is  dubbed  the  global  coupling  strength  between  the 
adjunct  and  the  force-driven  master  dynamic  systems.  The  loss 
factor  7;(®)  is  dubbed  the  induced  loss  factor;  the  loss  factor 
that  is  induced  in  the  master  dynamic  system  by  virtue  of  its 
coupling  to  the  adjunct  dynamic  system  [4-13,  15-17] .  In  the 
absence  of  coupling  fjM  is  identically  equal  to  zero.  In 

addition  to  the  virtual  loss  factor  defined  and  stated  in 

Eq.  (6),  an  effective  loss  factor  T/cW  combined  -  -  master 

+  adjunct  -  -  dynamic  systems  may  be  introduced.  This  effective 
loss  factor  is  defined  in  the  form 

Tie  (®)  =  ^  (^)]  '  7c  (^)  =  ®  ^  (®)]  '  W  +  (®)  ’ 

[n,(e»)/fi>£,(c>)]  =  7vH  =  7c(^)[l  +  ^^H]'  ^[{co)  =  [E,{m)lE,{a>)]  .  (7a) 

From  Eq.  •{7a)  the  ratio  between  the  virtual  loss  factor  and 

the  effective  loss  factor  TJg{a))  is  simply 

[;7»/7»]=[i+5»]>i  ,(7b) 
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where  the  ratio  {(o)  is  explicitly  expressed  in  Eq.  (6b).  It 
transpires  that  the  external  input  power  into  a  dynamic 

system  may  be  cast  in  the  form 

U^io))=S^(o))G{o))  !  g{o))={^/2)[v(co)/m]  ;  S^(o))=Acos  (oj) 

where  s^{o})  is  the  power  spectral  density  of  the  external  force 

drive,  V{o})  is  the  perceived  modal  density  and  (M)  is  the 
perceived  mass  of  the  dynamic  system.  Finally,  (Aco)  is  a  suitable 
frequency  bandwidth  centered  about  a  suitable  frequency  {p)^j  [3]  . 

The  perceived  quantities  are  those  perceived  by  the  force-drive. 

From  Eq.  (8)  the  external  input  power  Tlg{(o) ,  into  the  master 
dynamic  system  in  the  absence  of  coupling  to  the  adjunct  dynamic 
system,  may  then  be  estimated  in  the  form 

He (6))  S  ^ 

where  Vo{o})  and  {M^)  are  the  modal  density  and  mass  of  the  master 
dynamic  system,  respectively.  The  modal  density  V  (®)  and  the 

mass  {Mg,),  perceived  by  the  force  drive  in  the  absence  of 
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coupling,  may  be  modified  by  the  presence  of  the  coupling.  The 
modified  modal  density  and  the  mass  are  designated  by  and 

M^(ci)J ,  respectively.  It  is  to  be  assumed  that  the  modification 

is  related  to  the  averaged  modal  coupling  strength  in  the 

form 


V^W  =  [V^W+^^(a))V^(fi>)]  =  V^(6>)[l  +  3*W]  ,  (10a) 

Ml{(o)  =  M^[l  +  F{(MJMjgl{o>)]]  =  Mjl  +  F{{m^/mJ3l{o})}]  ,  (10b) 

where  \injan6.\in^)  are  the  modal  masses  in  the  master  dynamic 

system  and  in  the  adjunct  dynamic  system,  respectively,  and  it  is 
defined  that 


=  [v^{a))lV^{<o)] 

[mjmJ  [mjmJ 


, (lla) 


. (11b) 


Clearly,  Eq.  (10a)  is  less  problematic  than  Eq.  (10b)  [3,11]. 

Indeed,  the  functional  form  of  (f)  is  not  readily  definable  unless 
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is  small  compared  with  unity;  in  this  event, 
is  comparably  small  [11] . 

If  the  power  spectral  density  s ^{co)  of  the  external  force- 

drive  remains  intact  when  coupling  is  introduced,  then  from  Eqs. 
(8) -(10),  one  derives  for  the  external  input  power  11^ (<y)  the  form 

n^(fi;)s5^((i>)(;r/2)[v^W/Mj[l  +  3^(ty)][l  +  F{  / 772^)3^ (cy)}]'’  ; 

S ^{o})=  d^cos ^{co)  .  (9b) 

From  Eqs.  (3) -(10)  one,  finally,  obtains 

4  W = n:{fi>)4^  W = [i + 3*  W]  {[i + [t]^  (m)it]^  W)  3’  W]  [i + F  {(wp  (®)}]}'' 

,  (12a. 1) 

;^(®)  =  [l  +  3^(fi;)][l  +  F{(772^/m^)3^(®)}  ]"’  ,  (12b.l) 


(<y)  =  [  lo  (")/  Iv  («)]  =  [  1  +  {^5  lo  (<«)]  '  <  1 


. (12c) 
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Except  for  situations  of  special  design,  in  the  more  usual  noise 
control  designs,  the  modal  mass  ratio  is  small  enough  to 

render  small  compared  with  unity  even  when  the 

stored  energy  ratio  may  far  exceed  unity.  The  exceptions 

are  found  in  the  design  of  a  light  panel  (the  master  dynamic 
system)  that  is  destined  to  support  heavy  electronic  components 
(constituting  collectively  the  adjunct  dynamic  system)  in 
sections  of  space  crafts  [18].  By  and  large,  in  subsequent 

considerations  it  is  to  be  assumed  that  [(mjm„)%{co)\  is  small 
compared  with  unity  and,  therefore,  in  subsequent  considerations 

the  factor  +  is  approximated  as  equal  to  unity. 

In  this  vein,  Eqs.  (12a.  1)  and  (12b.  1)  are  approximated  in  the 
forms 

^W=[nlffl)/n°(<i>)]=[7,(‘»V7,W]=[i+3' w]>]  ,  (i2b.2) 


respectively,  where  the  effective  loss  factor  T]g{co)  is  defined  in 
Eq.  (7) .  Under  this  imposition  Eq.  (12)  already  poses  a  question 
of  significance:  May  an  adjunct  dynamic  system,  that  is  destined 
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to  be  passively  coupled  to  an  externally  force-driven  master 
dynamic  system,  be  appropriately  designed  to  achieve  a  credible 
noise  control  in  that  master  dynamic  system?  Three'  cases  are 
detailed  in  exemplifying  three  specific  but  diverse  answers  to 
this  question: 

Case  1.  In  this  case  one  assumes  a  priori  that  the  adjunct 
dynamic  system  is  merely  a  sink:  A  sink  is  a  dynamic  system  that 
absorbs  power  but  does  not  store  energy.  Thus,  in  this  case, 

3o(<y)=?0,  From  Eg.  (12b),  therefore,  7r{o))  is  identically  equal 
to  unity 

7r{co)  =  \  .  (13a) 

It  follows  that  when  the  adjunct  dynamic  system  is  a  sink  and  it 
is  attached  to  the  master  dynamic  system,  the  external  power 
injection  remains  unchanged  by  this  attachment.  On  the  other 
hand  the  loss  factor  (%)  that  characterizes  the  adjunct  dynamic 
system,  which  is  in  this  case  a  sink,  yields,  from  Eg.  (12),  the 
induced  noise  control  parameter  to  be 
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Here  the  induced  loss  factor  ?7/W  is  the  loss  factor  contributed 
to  the  master  dynamic  system  by  an  adjunct  dynamic  system  that  is 
a  priori  a  sink.  In  this  case  nM  is  commensurate  with  the  loss 
factor  that  is  also  assumed  a  priori  to  be  contributed  by  an 

attachment  to  a  sink.  From  Eq.  (13)  one  needs  recognize  that  if 
77/(<»)  is  to  be  finite,  cannot  be  selected  arbitrarily  small 

compared  with  unity. 

In  the  remaining  two  cases;  Case  2  and  Case  3,  the  adjunct 
dynamic  system  is  not  a  sink.  Indeed,  in  both  cases  it  is 

assumed  that  the  global  coupling  strength  rendered, 

by  design,  high  compared  with  unity;  3^(ft;)»l.  (Nonetheless,  as 

tacitly  assumed,  even  though  'v5„(®)»l,[(/Wj/Wj,)'Jo(6))]«l.  ) 

Case  2.  If  in  addition  to  rendering  -«5o(iy)»l,  the  loss  factor 
7]g{a>)  in  the  adjunct  dynamic  system  is  designed  to  highly  exceed 
the  loss  factor  that  is  inherent  to  the  master  dynamic 

system;  then  from  Eq.  (12)  one  obtains 
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which  describes  a  beneficial  noise  control.  It  needs  to  be  said 
in  this  connection  that  were  the  ratio  7l{p^  of  the  external  input 
power  assumed  to  be  equal  to  unity,  as  stated  in  Eq,  (5a)  ,  the 
apparent  noise  control  achieved  under  this  (false)  assumption, 
would  be  even  more  beneficial  than  that  estimated  in  Eq.  (14a)  ; 
i .  e . , 

!(<»)=>  4 (®)£[i  +  {;7,(<B)/7„(43;(ffl)r'<  [riA<»)iriA‘4  «  i  • 

Case  3.  If  in  addition  to  rendering  ^  \q))»1,  the  global 

coupling  strength  can  be  rendered  high  enough,  such  that 

even  if,  by  design,  is  still  in 

excess  of  unity,  then 

4(co) s[riAcAiri, (<»)]>  ^  -dSa) 

Equation  (15a)  describes  a  noise  control  reversal.  Again,  were 
the  ratio  7r{o})  of  the  external  input  power  assumed  to  be  equal  to 
unity,  as  stated  in  Eq.  (5a),  the  noise  control  reversal  would 
not  emerge;  namely,  under  this  (false)  assumption 


13 


,  (15b) 


which  is  3  bensficial  noise  control^  thus,  contnaciicting  the 
estimated  noise  control  reversal  quoted  in  Eq.  (15a) .  (Again  it 

is  recalled  that  in  this  consideration  is  assumed 

to  be  small  compared  with  unity.) 

A  corollary  to  cases  2  and  3  follows.  Were  the  adjunct 
dynamic  system  a  loss  factor— wise  similar  to  the  master  dynamic 
system,  in  the  sense  that  will  find  that  Eqs . 

(12a)  -  (12c)  assume  the  forms 

7t{a))  =>  [l  +  3o(o)]  >  1 

<?oW=[i+^oWr*  <1 

respectively.  Thus,  when  the  coupling  introduces  no 

noise  control  benefit;  the  coupling  is  neutral. 


, (16a) 


,  (16b) 
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In  addition  to  the  obvious  roles  played  by  the  indigenous  loss 
factors  rjoio))  and  77,  (ft?)  (and  the  modal  masses  and  ,  in 

the  master  and  in  the  adjunct  dynamic  systems,  respectively,  in 
the  determination  of  the  noise  control  parameter  <^{oy) ,  the 

crucial  role  played  by  the  global  coupling  strength  is 

clear  in  Eg.  (12),  especially  when  the  adjunct  dynamic  system  is 
not  merely  a  sink.  One  recalls  that  in  a  sink 

and  ?J^{a))^l{co)^?Jj{co) .  The  nature  and  the  composition 
of  a  non-zero  3o(<y)  is,  therefore,  investigated  next. 
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III.  Global  and  modal  coupling  strengths  and  masses 

As  already  intimated  in  Egs.  (10)  and  (11),  the  modal  coupling 
strength  is  related  to  the  global  coupling  strength 

in  the  form 


E^{o))=AcoVM^o{^)  '  ^,W  =  A^»V,(o)^,W  ; 

(0)  =  Afl)V„  (o)  =  (Ao / 0)^  )[(OoV„  (fi>)]  ;  (0)  =  ^(oV^  {o>)  =  {Aco /  o)^  )[<a„V^  (a))] 

,  (18a) 

where  V^(®)and  S„{a})  and  V^(i»)and  are  the  modal  density  and 

the  modal  stored  energy  in  the  master  and  in  the  adjunct  dynamic 
systems,  respectively,  (A<y)  is  a  suitable  freguency  bandwidth 
centered  about  (ft>^)[l-3].  Averaging  a  la  Skudrzyk  is  implied  in 

the  likes  of  Eqs.  (17)  and  (18)  and  in  some  subsequent  equations 
[14]  .  [cf.  Appendices  A  and  B.j  In  order  to  proceed  further,  it 
is  necessary  to  define  more  precisely  the  modal  density  V^{co)  of 

the  adjunct  dynamic  system.  For  this  purpose  the  distribution 
X /,X  ^ !  (O^ ,  is  required,  where  (<»,.)  is  the  modal  resonance 
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frequency  of  the  (r)th  mode;  the  modes  are  assumed  to  be 
sequentially  indexed.  If  the  index  (r)  is  given  a  continuous 
connotation  then  the  continuity  and  sequentiality  of  the 
distribution  may  be  defined  in  the  form 

X{r)=\_(o{r)l  co^  ;  x{r  +  s)>  x{r  -  e)  ;  £>0  .  (18b) 

One  may  then  define  the  corresponding,  local  modal  density  in  the 
form 


kv(r)]  =  [ay(r)/ar] 


.  (18c) 


With  this  definition  in  place,  mean-value  estimates  of  the 
induced  loss  factor  JJjioj)  may  be  determined  [3,  13,  15,  16]  .  The 

result  is 


T], {a,)s(nl2)(ala,J[a„V,(co)](m,  / WJcW 


.  (19a) 


The  coupling  factor  C{co) ,  in  Eq.  (19a),  is  expressed  in  terms  of 
the  normalized  coupling  coefficients  in  the  form 
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c(<»)  =  [l  +  m^(fl;)]  '  c„(<y)  ;  C^(<y)  =  [  {ni^W+ac(®)}%{^(®)}^]  ,  (20a) 


m^W=[m^W/mJ 


a ^  (o))  =  [k^  {co)/{(D^m^  )  ]  ;  g{o))  =  [g  (o))/  to  W  J  .  ( 2  Ob ) 


In  Eq.  (20)  the  vector  ^W}  defines  the  mass,  the 

stiffness  and  the  gyroscopic  coupling  coefficients,  respectively. 
Conveniently,  one  may  categorize  the  coupling  factor  C{o))  from- 
strong-to-weak  in  the  form  [13,  15] 


(l/;r);  strong  couplings 

, (21a) 

C{a))= 

{tcXO^);  moderate  couplings 

,  (21b) 

(lO"^);  weak  couplings 

. (21c) 

One 

may  utilize  the 

relationships  between  the 

global  masses 

(m) 

and  {.M ^  to  the 

corresponding  modal  masses 

{mj  and  {mj; 

namely 


[N^{a>)mjN^{(D)mJ  =  {MJM„)  ;  M^=N^{co)m„  ;  =N^{(o)m^  ,  (22) 


and  Eq.  (18a)  to  recast  Eq.  (19a)  in  the  form 
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r],(o,)  =  (nl  2)(ffl  /«„)’[<»„  V„  /  A/ Jc(®) 


.  (13b) 


Again,  it  is  emphasized  that  Tlj{o))  is  the  induced  loss  factor  that 
defines  the  net  power  flow  ns(fi>)  from  the  master  dynamic  system 
to  the  adjunct  dynamic  system  in  terms  of  the  stored  energy 
in  the  former;  i,e.,  Eq.  (6a),  in  part  states 

n,(6>)  =  /77W[®^oW]  -(23) 

Now  that  the  induced  loss  factor  T] ^{co)  has  been  defined  and 
stated  in  some  detail,  attention  is  to  be  focused  on  the 
definition  and  statement  of  the  modal  coupling  strength  ?»■ 

From  Eqs.  (6b)  and  (17)  one  finds  that  to  derive  off 

one  needs  define  the  indigenous  loss  factor  in  the  adjunct 

dynamic  system  [13,’  15,  16]  .  For  this  purpose  it  is  convenient 

to  express  in  terms  of  the  modal  overlap  parameter  h^{o)) 

[13,  15,16].  The  expression  is 

rj^  {(o)  =  {o})[co  {o))]  ^  {co)[o)^V^  {a)]  .(24) 
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where  the  modal  density  V^(o)  conforms  to  its  definition  in  Eq. 
(18c)  .  The  convenience  stems  from  the  role  that  is  played  by 
b^[(o)  in  relating  frequency  widths  of  modes  to  the  frequency 

separation  between  adjacent  modes  [3].  In  that  sense  b^{a)^ 
determines  the  degree  of  damping  that  represents;  namely 

«  1,  light  damping 

=  1,  moderate  damping  .  (25) 

»  1,  heavy  damping 

Indeed,  although  the  values  of  stated  in  Eq.  (19)  appear  to 

be  independent  of  the  modal  overlap  parameter,  this  independence 
is  not  free.  The  derivation  of  T]j{(o)  in  the  form  stated  in  Eq. 

(19)  is  predicated  on  a  priori  assigning  a  continuous  connotation 
to  the  indices  that  identify  the  modes  in  the  adjunct  and/or  the 
master  dynamic  systems.  The  assignment  allows  the  summations 
over  modes  in  the  evaluations  of  the  induce  loss  factor  to  be 
replaced  by  integrations  over  the  continuous  indices.  It 
transpires  that  this  process  yields  average  (a  la  Skudrzyk)  values 

for  TIj{g>)  when  b^{a))  is  less  than  unity.  (The  summations  in  this 

instance;  i.e.,  when  b^{ci})  is  less  than  unity,  yield  induce  loss 
factors  that  undulate  as  functions  of  the  normalized  frequency 
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{col co^)  .  The  excusions  of  the  undulations  from  the  mean-values 

are  the  more  pronounced  the  lower  are  the  values  of  b^{co) . 

[Appendix  A  and  B.])  On  the  other  hand,  when  b^{co^  exceeds  unity, 
the  evaluation  of  the  induced  loss  factors  by  summations  over  the 
individual  modes  substantially  match  those  values  derived  by 

integrations  over  the  continuous  indices.  In  this  sense  bg{coi) 
plays  a  pivotal  role  in  the  interpretation  of  the  data  yielded  in 
the  subsequent  developments  of  this  thesis.  From  Eqs.  (6b), 

(17),  (19)  and  (24)  one  obtains 


gl{a)sU/2){o,/o,J  [<i)„  V.  /  mJlc{o)/b,  (»)] 


, (26a) 


where 


[n(®)/V„W]  =  [w,(®)/Af„W]  .  (27) 

The  induced  loss  factor  stated  in  Eq.  (19)  ,  and  the 

modal  coupling  strength  <^^{co),  stated  in  Eq.  (26),  are  exemplified 
for  the  complex  dynamic  system  sketched  in  Fig.  3.  This  complex 
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dynamic  system  comprises  a  simple  harmonic  oscillator  for  the 
master  dynamic  system  and  a  number  of  simple  harmonic  oscillators 
that  collectively  constitute  the  adjunct  dynamic  system.  The 
coupling  elements  are  uniformly  composed  of  individual  or 
combined  mass,  stiffness  and  gyroscopic  components  [13,  15,  16]  . 
Only  the  master  dynamic  system  is  externally  force-driven;  the 
satellite  oscillators,  that  constitute  the  adjunct  dynamic 
system,  are  neither  coupled  to  each  other  nor  externally  force- 
driven.  The  so  described  complex  dynamic  system  is  amendable  for 
precise  definition  in  terms  of  the  individual  oscillators, 
including  the  master  oscillator,  and  the  individual  couplings  of 
the  satellite  oscillators  to  the  master  oscillator.  [Appendix 
B.]  Moreover,  the  so  prescribed  complex  dynamic  system 
facilitates  the  computations  of  the  designated  quantities,  'IJjico) 

and  rendering  the  results  of  these  computations  easy  to 
interpret . 
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IV.  Computations  of  the  induced  loss  factor  Tjj{o))  and  the  modal 
coupling  strength 

To  facilitate  the  computations  a  more  definitive  statement 
must  be  introduced  regarding  the  complex  dynamic  system.  Thus, 
from  Fig.  3  one  surmises  that  is  equal  to  unity  and, 

therefore,  mg=Mg.  [cf.  Equation  (22).]  In  addition,  the 
resonance  frequency  of  the  master  dynamic  system  in  isolation  is 
designated  {cOg).  The  harmonic  oscillators  pertaining  to  the 
adjunct  dynamic  system  are  distinguished,  in  this  case,  as 
satellite  oscillators.  The  number  Ng{o})  of  satellite  oscillators 
may  exceed  unity;  iV^(®)>l.  In  addition,  the  resonance  frequency 
of  the  (r)th  satellite  oscillator  is  designated  (^r)  ■  '^he 
normalized  resonance  frequency  distribution  Jf(r)  of  the  satellite 
oscillators  may  be  defined  in  the  form 

!  (Og)=\co{r)l  o)g\=  X{r)  ;  \<r<N^  ,(28) 

where  (r)  indexes  the  (r)th  satellite  oscillator.  [cf.  Eq. 
(18b).]  The  functional  form  of  X{r) ,  as  a  function  of  (r),  may  be 
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given  a  continuous  connotation;  e.g.,  as  exemplified  in  Fig.  4a. 
In  this  figure 

+  +  ;  r  =  ^-^  '  ,  (29a) 

and  Nj=21  [15,16].  Both,  the  discrete  and  the  continuous,  forms 

of  X(r)  are  depicted  in  Fig.  4a.  [cf.  Appendix  A.]  Also,  it  is 
imposed,  in  this  example,  that  as  many  satellite  oscillators  are 
with  resonance  frequencies  that  are  less  than  the  resonance 
frequency  {o)o),  as  with  those  that  exceed  {(Og)  [13]  .  Indeed,  when 
the  number  of  satellite  oscillators  is  odd,  the  resonance 
frequency  of  one  of  the  satellite  oscillators  matches  that  of  the 
master  oscillator;  for  that  oscillator,  by  definition, 

+  =(l/2)  and  ^(r)  =  l  [13].  Similarly,  the  loss  factor  , 

that  is  associated  with  the  (r)th  satellite  oscillator,  may  be 
cast  in  the  form 


rj^^ri(r)=b(rtx(rj\dX(r)ldr]^b(rtrnN,)[x(rf  ,  (29b) 

where  b{r)  is  a  designated  localized  modal  overlap  parameter; 
localized  at  and  in  the  vicinity  of  the  resonance  frequency  (cOj) 
of  the  (r)th  satellite  oscillator  [13].  Employing  Eq.  (29a), 
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with  N^=27  in  Eq.  (29b),  J]{r)  is  exemplified  in  Fig.  4b;  in  this 

example  three  constant  values  of  b{r)  are  used;  namely,  Z)(r)=(0.1)  , 
(2.0)  and  (10)  [15,  16].  Again,  both  the  discrete  and  the 
continuous  forms  of  7j{r)  are  depicted  in  Fig.  4b.  [cf.  Appendix 
A.] 


Using  the  expression  for  the  normalized  resonance  frequency 
distribution  X{r) ,  as  stated  in  Eq.  (29a),  one  derives  from  Eq. 
(18c),  the  expressions 


[<y„V^(^y)]  =  N^ia)  (2lr)(o)o  Icof 


,  (30a) 


{/^(Olco^)  '  =(2lr){(olco^Y  =[(o^V^{o))\ 


-1 


, (30b) 


where,  again,  AC(<y)  =  l.  From  Eqs.  (19)  and  (30)  one  obtains 


rij{a))  =  (7rl2)(2/r)(M^  I M )c{(o) 


(31) 


The  values  for  the  induced  loss  factor  ,  as  a  function  of  the 

normalized  frequency  {colcoY),  in  the  appropriate  range  of 
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frequency,  are  depicted  in  Fig.  5. 
frequency  is  given  by 


The  normalized  range  of 


[AW/tt>Js[l-(r/2)]  '^^-[l  +  (r/2)r^^  =(r/2)  .  (32) 

Again,  in  Fig.  5,  iVj=27and  the  global  mass  ratio  is  set 

ecjual  to  on©  t©nth  (1/10)  •  Thsr©  ©re  thro©  curves  in  Fig.  th© 

solid  curve  corresponds  to  a  strong  coupling  with  C(<o)=l,  the 

dash  curve  corresponds  to  a  moderate  coupling  with  C(fi?)=3xl0  and 

the  dash-dot  curve  corresponds  to  a  weak  coupling  with  C(to)  =  10  . 

The  normalized  range  of  frequency  in  Fig.  5  is  set  by  }'  =  0.6. 
Clearly,  the  induced  loss  factor  increases  with  an  increase  in 

the  global  mass  ratio  with  an  increase  in  the  coupling 

factor  C{o))  and  with  a  decrease  in  the  normalized  range  of 

frequency;  namely,  a  decrease  in  (^);  notwithstanding  that  (^)  must 
be  chosen  less  than  unity.  It  may  be  of  interest  to  contrast  the 
induced  loss  factor  TJ^ico)  with  the  corresponding  assigned  values 

for  the  indigenous  loss  factor  TJ^{co)  in  the  adjunct  dynamic 

system.  The  expression  for  is  derived  from  Eqs.  (24)  and 

(30)  to  be 
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Vs  (®) = ^s  (®)(® '  f  Ir '  '^^s  (®)] 


.  (33) 


The  loss  factor  V depicted,  as  a  function  of  {colco^),  for 
7  =  0.6  and  N^{co)  =  21 ,  in  Fig.  6.  The  three  curves  in  Fig.  6  pertain 

to  Z>^(g;)  =  0.1  (solid),  b^{a))=2Xi  (dash)  and  b^=\0  (dash-dot).  The 
results  in  Figs.  5  and  6  may  be  compared  as  a  bench  mark. 

Using  Eqs.  (30)  and  (33)  in  Eg.  (26)  one  obtains  for  the  modal 
coupling  strength  the  equivalent  expressions 

f'(<i))  =  [;r/2Z>,(<i.)](<i>„/<»f(2/,)^K(ffl)](OT,/W„)[c(®)]  ,  (34a) 

$-'(®)s(>r/2)U„/<»)'(2/,f(w,/V,)[cW/4(<i>)]  ;  N„(a)=\  .  (34b) 


The  values 

for  the  induced  loss  factor. 

as  a  function  of 

the 

normalized 

frequency 

{cd/o)^),  in  the 

appropriate 

range 

of 

frequency. 

as  defined 

in  Eq.  (32) ,  are 

depicted 

in  Fig. 

7. 

Again,  in  Fig.  7  the  global  mass  ratio  is  set  equal  to  one-tenth 
(1/10);  (A/'^/Mp)  =  0.1 ,  and  the  normalized  range  of  frequency  is  set 
by  7  =  0.6.  Clearly,  as  is  the  induced  loss  factor 
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modal  coupling  strength  increases  with  increase  in  the 

global  mass  ratio  and  the  coupling  factor  C{o)) .  Also,  as 

does  TJjifi)),  increases  as  (r)  decreases.  Significantly, 

however,  the  modal  coupling  strength  (ft))  is  inversely 

proportional  to  the  modal  overlap  parameter  .  As  Eg.  (34b) 

indicates  one  may  define  [c(<y)/Z?j(o)]  as  the  modified  coupling 

1 ,  . 

factor;  the  modified  coupling  factor  C  (co)  exceeds  the  coupling 
factor  if  is  less  than  unity  and  is  less  than  the  coupling 

factor  if  b^ito)  exceeds  unity.  In  each  sub-figure  of  Fig.  7  three 
curves  are  depicted;  the  solid  curve  pertains  to  a  b^io^  that  is 
equal  to  one-tenth  (1/10),  the  dash  curve  pertains  to  a  b^{G))  that 

is  equal  to  two  (2)  and  the  dash-dot  pertains  to  a  b^{a))  that  is 
equal  to  ten  (10) .  Figure  7a  depicts  strong  coupling  with  a 
coupling  factor  C(o)  that  is  equal  to  unity;  C(©)  =  1,  Fig.  7b 

..  -2 

depicts  moderate  coupling  with  C{ct>)  =  3x\0  and  Fig.  7c  depicts 

weak  coupling  with  C(e>)  =  10^.  A  major  observation  emerges  from 
Figs.  7a-c;  the  observation  is  discussed  briefly  in  the  next 
section. 
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V.  Validity  of  statistical  energy  analysis 

An  observation  emerges  when  the  modal  coupling  strength 

sea  ,  V 

^  {(o)  is  estimated  via  the  statistical  energy  analysis  (SEA) 

rather  than  estimated  in  the  manner  described  herein;  e.g.,  as 
stated  in  Eg.  (34)  [13]  .  In  terms  of  SEA  the  net  power  11^ (^y) 

that  invades  the  adjunct  dynamic  system  is  given  by 

n J  W  =  t]„  (o>XaE^  (»)]  -  ri^  (®)[®£,  (®)]  ; 

where  77  ((o)  and  T?  (co)  are  the  coupling  loss  factors  from  the 
•so^  '  ' os^  ' 

master  dynamic  system  to  the  adjunct  dynamic  system  and  vice 
versa,  respectively  [3] .  From  Eq.  (35)  one  obtains 

?rw=7„(®H(®)+'7„W]''  -  (36«) 

or  equivalently 
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[  n,  (®)/  (0)]"“  =  {®) = %  (®)[%  H + '7..  (®)]'' = 3r  (<»)  '• 

3r(®)=[Af.w/-w„w]$'rw  ' 

where  the  superscript  (sea)  indicates  that  the  estimates  so 
superscripted  and  those  enclosed  in  the  square  brackets  are  a  la 
SEA  [3]  .  It  is  clearly  evident  from  Eq.  (36)  that  the  modal 

coupling  strength  {(o)  must,  by  definition,  lie  below  unity; 

namely 


1  •  (37) 

Equation  (37)  is  a  tenet  of  the  statistical  energy  analysis  (SEA) 
[3,  19]  . 


For  the  complex  dynamic  system  depicted  in  Fig.  3  and  in  terms 
of  the  analysis  presented  herein  this  restriction  on  the  modal 
coupling  strength  may  be  stated,  from  Eqs.  (34)  and  (39),  to  be 


bs(a>)>[b^{a>)]M  /of  {4/ /)(m^  IM^)c{g)) 


,  (38) 
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where  [bA-X  designates  the  minimum  value  of  the  modal  overlap 
parameter  that  is  needed  to  validate  the  statistical  energy 
analysis.  In  Fig.  8  the  values  of  [Z)^(<y)]^,  as  a  function  of 
{colcOg),  is  depicted  for  three  values  of  the  coupling  factor 

C(<y);  C(f»)  =  l  (dash- dot ),  C(<y)  =  3x10  (dash)  and  C(6))=  10  (solid).  It 
is  observed,  in  Fig.  8,  that  (SEA)  is  validated  for  values  of 
b^{G))  that  are  less  than  unity  only  for  coupling  factors  that  are 

small  compared  with  unity.  It  is,  thus,  speculated  that  the  use 
of  (SEA)  may  be  strained  when  dealing  with  strong  and  even  with 
moderate  couplings.  The  validation  of  (SEA)  may  be  called  to 
question  when  the  loss  factors  that  are  associated  with  the 
adjunct  dynamic  systems  assume  arbitrarily  low  values  and  the 
coupling  factors  (and  the  global  mass  ratios)  assume  high  values. 
In  this  connection  one  is  reminded  that  undulations  in  the  non¬ 
average  values  in  the  response  quantities  of  a  complex  dynamic 

system  are  suppressed  only  when  b^{a}^  exceeds  unity,  [cf.  Appendix 
B.]  However,  the  absence  or  the  presence  of  undulations  in  that 


non-average 

response 

do 

not 

bear  on  the 

criterion  for  the 

validity  of 

(SEA) ; 

(SEA) 

is 

validated  as 

long  as  the  modal 

overlap  parameter  b  ((5>)  exceeds  the  corresponding  minimum  value  of 
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this  parameter;  •  Indeed,  the  employment  of  (SEA) 

may  be  valid  for  values  of  both  and  b^{(o)  that  may  be  less 
than  unity.  Equation  (38)  makes  clear  that  the  validity  of  (SEA) 
merely  requires  that  b^(o^  exceeds  •  The  validity  does  not 

demand  that  undulations  be  absent  in  the  non-average  values  of 
the  response  quantities.  [cf.  Appendix  B.]  Of  course  (SEA) 
yields  largely  average  response  quantities;  which  are 
commensurably  averaged  a  la  Skudrzyk  [3,  14] . 
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Appendix  A.  Discrete  and  Continuous  Distributions 


The  expression  for  the  resonance  frequency  distribution  •  is 
formally  stated.  In  the  master  dynamic  system  the  expression  is 

=  ,  (Ala) 

and  in  the  adjunct  dynamic  system  the  expression  is 

{(o^lco^)^[co{r)lco^]  =  xXr)  ,  (Alb) 

where  (n)  and  (r)  designate  the  (w)th  and  the  (r)th  modes  that 
reside  within  the  frequency  bandwidth  (Ag>)  centered  about  the 

frequency  (o^),  in  the  master  and  in  the  adjunct  dynamic  systems, 
respectively,  and  X^(n)  and  Xg(r)  are  the  corresponding  modal 

distributions  in  the  frequency  band  defined  by  {Aci?,£i}J  [3,  13,  15, 
16]  .  In  Eq.  (Al)  the  quantities  and  parameters  are  stated  in  a 
manner  that  allows  for  a  conceptual  smear  from  the  discrete-to- 

the  continuous  in  the  frequency  band  .  As  an  example,  the 

quantities  X^(n)  and  are  chosen  to  be  in  the  forms 

Jir„(«)=[l  +  {(l  +  W„)-2«}(r„/2W„)]-'  ;  ^,,=0.3  ;  W»>1  ,  (A2a> 
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(r)  =  [1  +  { (1  +  W,  )  -  2r}  O'  /2Ar. 


;  /  =  0.6  ;  Ng  (co)  >  1 


.  (A2b) 


so  that  they  reasonably  span  the  same  frequency  band  {^(0,(0^} 
[13,15].  [cf.  Fig.  4a.]  In  this  connection,  one  may  further 
define  from  Eqs.  (Ala)  and  (Alb)  the  loss  factor  distributions  in 
the  master  and  adjunct  dynamic  systems 


, (A3a) 


and 


Tlr-‘‘nA’-)=b,('-)UA'V  ^X,(r)iar\  ,(A3b) 

respectively,  where  b„(n)  and  b,(r)  are  the  respective  local  modal 
overlap  parameters.  From  Eqs.  (A2)  and  (A3)  one  obtains 

rj^  («)  =  {n)[x^  (n)]  IN^)  ;  (®)  =>  b„  {a))(a)  I  a)  ^)[y^lN^  (®)]  ,  ( A4  a ) 

;  riXco)^bX(3>)(o>l(Oof[rl2NX(o)\  ,  (A4b) 

where  the  conversion  from  the  index-domain  into  the  frequency- 
domain  is  governed  by 
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(cdIco^)^  X^{n)  ;  (<y/<yj  =>  X^(r)  ,  (A5) 

and  it  is  observed  that  the  frequency  bandwidths  are  defined  by 
the  values  of  X^{n)  and  X^(r)  at  the  end-points;  namely 

X„(l/2)  =  [l  +  (r„/2)r'  ;  2r„{Af„+(l/2)}  =  [l-(r„/2)r'  ,(A6a) 

Xll/2)  =  [l+(r/2)]‘‘'^  ;  X,{N,+(\l2))-=[\-(rl2)T''^  .  (A6b) 

[cf.  Fig.  4b.]  These  limits  are  commensurate  with  a  frequency 
bandwidth  A(<y)  in  the  form 

[A(ft))/<y^]  =  70=  (yl'2) ;  2yo=r  =  0-6  •  (A7 ) 

On  the  other  hand,  the  values  of  the  average  loss  factors 
{TJ^  and  are  estimated  from  Eqs.  (A3)  and  (A6)  to  be 

_1  No  _1  {A^o+(l/2)} 

»7<,=(JV<,)  Z'7<,W=>(A'.)  J  ; 

n  (1/2) 

no=<Vo(’‘)>=(Nj'  \b,(nlx,(nj'  dX,{’’)  =  K^ol  N,)  ,  (A8a) 
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r  (1/2) 

7,  =<%('-)>=(^.)'‘ J*.('-fe('')r'‘if,('')s4(r,^2^.)  • 


The  replacement  of  a  summation  by  an  integration,  as  is  stated 
cavalierly  in  Eq.  (8),  needs  to  be  properly  qualified.  The 
extrapolations  and  the  interpolations  of  the  summand  to  establish 
the  integrand,  as  the  indices  (r)  and/or  (n)  are  given  continuous 
forms,  must  render  the  latter  a  fairly  smooth  function  in  the 
continuous  variables  (r)  and/or  (n)  ,  respectively.  This 
smoothing  of  the  integrand  does  not  only  facilitate  the 
replacement  of  the  summations  by  integrations;  e.g.,  as 
exemplified  in  Eq.  (8),  but  this  smoothing  is  an  averaging 
process  similar  to  the  mean-value  averaging  proposed  by  Skudrzyk 
[14] .  To  illustrate  this  process  and  render  it  more  concrete  one 
may  proceed  to  analyze  Eq.  (A8b)  in  more  detail.  Concentrating 
on  the  first  of  this  equation,  one  may,  by  extrapolations  and 
interpolations,  insert  a  step  in  the  form 


r=l 


Y  T]Xr  +  {r'/K)] 

r'  =  (-iV;/2) 


(A9a) 
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where  (r')  is  a  running  index  and  N'^  is  as  large  compared  with 
unity  as  is  suitable  to  render  the  summand,  in  this  equation, 
smooth  as  a  function  of  [r  +  (r7iV')]  with  (r)  given  one  of  its 
sequential  integer  values.  Then,  in  the  limit,  if  the  summations 
in  Eq.  (A9a)  are  absolutely  convergent  (not  conditionally 
convergent)  one  may  state  the  equivalence 


iNjVfi  E 

r=l  r'=iN's/2) 

^  (A'.)"' 

where  the  index  (r)  in  the  integral  is  piece-wise  continuous  with 
a  grain  that  is  determined  by  the  largeness  of  (iV')  compared  with 

unity.  Clearly,  iK)  may  be  set  as  high  as  smoothing  may  require 
and  beyond.  In  the  limit  the  summand,  in  the  sum  over  (r)  from 
r=l  to  ,  is  smooth  and  the  summation  qualifies  for 

replacement  by  integration.  Notwithstanding  that  were  the 
indices  (r)  and/or  (n)  render  a  priori  continuous  in  the  summand, 
would  have  resulted  in  an  integrand  that  is  smooth,  so  that  the 
replacement  of  the  summation  by  an  integration  would  be  naturally 
qualified  [6].  It  is  provided  that  within  the  limits: 

\<r<Ng  a.ndil<n<N^  the  modal  overlap  parameters  b^{r)  and  b^{n) 
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and  the  resonance  frequency  distributions  X^{r)  and  X^{n)  are 

reasonably  smooth  functions  of  the  continuous  (r)  and  (n) , 
respectively;  e.g.,  as  represented  in  Fig.  4.  The  replacing  of 
the  summations  by  the  integrations  in  Eg.  (8)  is  then  naturally 
qualified. 

Attention  is  now  focused  on  the  induced  loss  factor.  It  is 
further  developed  that  the  discrete  expression  for  the  induced 

loss  factor  TJj(n)  and  its  continuous  counterpart  7Jj{n) ,  both 
pertaining  to  the  (n)th  mode  in  the  master  dynamic  system,  may  be 
cast,  respectively,  in  the  forms 


N, 


%  fJiM !  %W=W)  Z 


iK/2) 

E  ^/[  n,r  +  {r'/N')]  ,  (AlOa) 

r'=(-JV;/2) 


and  their  averages,  over  the  modes  in  the  master  dynamic  system, 
are  given  by 


Tjf  =<t}f(n)MN,)-'Z  nf(n)  ■ 

n 

No  iNo/2) 

r],=(Ny(Ky"L  z  .  (Alla) 

n  n'={-N'ol2) 
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[cf.  Eq.  (A9).]  If  (a^')  and  (a^^)  are  assigned  large  enough  values 

and  the  extrapolations  and  interpolations  are  properly  chosen  to 
render  the  suitimand  smooth,  the  summations  may  qualify  to  be 
replaced  by  integrations.  Then,  the  values  of  the  average 
induced  loss  factors  T]j{n)  and  TJj  are  estimated,  from  Eqs .  (AlOa) 

and  (Alla),  to  be 

{TV, +(1/2)} 

rii{n)^  \rij{n,r)dr  ,  (AlOb) 

(1/2) 

{iV„+(l/2)}{{/V^+(l/2)} 

TIj=>{N^)  j  J  ?Jj{n,r)drdn  .  (Allb) 

(1/2)  (1/2) 

Again,  the  replacement  of  a  summation  by  an  integration,  as 
stated  in  Eqs.  (AlO)  and  (All),  is  thus  procedurally  qualified, 
[cf .  Eq.  (A9) . ] 

Once  again,  it  may  be  advantageous  to  institute  a 
transformation  of  variables  in  Eqs.  (AlOb)  and  (Allb);  namely 

%(")=■  /%('<.'•)[/.  -  (MOc) 
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where  the  limits  on  the  integrals  over  Xg{r)  and/or  over  Xg{n)  are 
conveniently  omitted,  but  are  obviously  those  expressed  in  Eg. 
(A6)  and 

;  f,{rh{dXAr)ldr]-'  •  (A12) 

Thus,  the  utilization  of  extrapolations  and  interpolations  in  the 
procedure  to  replace  the  summations,  in  Eqs.  (AlOa)  and  (Alla), 
by  the  integrations,  in  Eqs.  (AlOc)  and  (Allc) ,  entails  rendering 
the  integrands  in  these  integrations  fairly  smooth  functions  of 
the  continuous  forms  of  and/or  respectively.  [cf. 

Fig.  4a  and  Eq.  (A9b) . ]  The  assurance  that  this  smoothness  can  be 
implemented,  with  respect  to  and  then  readily  made. 

For  example,  from  Eqs.  (A2)  and  (A12)  one  obtains 

/,(«)= ,(A13a) 
/,(r)=(2Af,/r){jf,(r))-’  > 
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which  are  reasonably  smooth  within  the  chosen  frequency  bandwidth 
indicated  in  Eqs .  (A6)  and  (A7).  It  is  to  be  understood  that,  in 
Eqs.  (AlOc)  and  (Allc) ,  where  appropriate,  the  vector-variable 
{n,r]  needs  to  be  replaced  by  .  Here  the  assumption  is 
that  these  two  two-vectors  are  properly  related,  whether  they  lie 
in  the  discrete  domain  or  in  the  continuous  domain.  [cf.  Eq. 
(A5).]  Unfortunately,  without  further  investigation  no  assurance 
of  smoothness  can  be  readily  given  with  respect  to  the  integrand 
factor  Tlj{n,r)  in  Eqs.  (AlOc)  and  (Allc).  The  smoothness  or  lack 

of  smoothness  in  7/ ("’'■)  is  examined  in  the  next  appendix;  i.e.,  in 
Appendix  B. 
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Appendix  B  Smoothing  of  the  Induced  Loss  Factor  Tl,(n.r) 

The  equation  of  motion  that  governs  the  induced  modal  loss 
factor  rij(n,r),  between  the  («)th  mode  in  the  master  dynamic  system 
and  the  (/‘)th  mode  in  the  adjunct  dynamic  system,  may  be  expressed 
in  the  form 


77,  («. r)  =  -ImM  (m,  /»!„)[ {1  -  (z. )' (1  +  i T/r )}  { -  (z„ (l  +  )}  -  (q„  I yf] 

[(i+m„)-(z^)"(i+i77„„)]''}  .  (B1 


where  a  similarity  condition  is  imposed  in  the  form 

(z,f  ^aAxAr)iyf  ■  {z,rf  =aJxAr)/y?  ;(z„)"  =(«,  ya^lx,(r)iyf  ■ 

y  =  (ci>lm,)  ;  (z^)^(l  +  '77„^)=(Z,)^(l  +  /?7r)+(Z»-f  +  '  (B2a) 

and  the  coupling  parameters  (m^)and(^^)  are 

m„=(m^im,)  ;  (q„f  =  •  (B2b) 

The  parameters  that  describe  the  adjunct  dynamic  system  -  -  mass- 
stif  fness-dashpot  -  -  are  (/W,.),  {Clf,. ,  Xj(r)}and  (T/,.)  and  those  that 

describe  the  coupling  between  the  master  and  the  adjunct  dynamic 
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systems  -  -  mass-stif fness-dashpot-gyroscopic  -  -  are 
(m„,),  {«„,,X,(r)},(7„,)and  (g„,).  The  similarity  imposed  on  the 

stiffness  associated  with  the  mass  iplr)>  both  with  respect  to  the 
(r)th  mode  and  with  respect  to  the  coupling  of  that  mode  to  the 
(?7)th  mode  in  the  master  dynamic  system,  allows  one  to  incorporate 
the  couplings  in  the  description  of  the  adjunct  dynamic  system. 
This  incorporation  is  a  useful  interpretative  ploy  [13,  15]  .  The 

resonant  nature  of  ?Jj(n,r) ,  stated  in  Eq.  (Bl)  is  obvious.  The 

modes  are  defined  in  the  -domain  by  requiring  the  quantity 

[(l  +  ]  to  vanish.  Then  the  width  of  the  mode  in  the 

-domain  is  determined  by  [0  + as  the  local 
modal  overlap  parameters  Q^nrr^  below  unity,  the  modes  are 

distinguishable  and  the  distribution  of  T]j{n,r) ,  for  discrete 
values  of  the  index  (r)  and  a  fixed  value  of  the  index  (n) ,  is 

not  smooth  as  a  function  of  (^y/^y^)  in  .  The  presence  of 

modes  in  rij(n.r),  as  a  function  of  {cotco^  in  {Aco/ CO^) ,  for  discrete 

(r)  and  (n) ,  is  readily  discernible.  Indeed,  the  presence  of 

modes  in  the  summand  TJj{n,r) ,  in  this  case,  is  even  discernible  in 

d,  V 

the  summ  stated  in  the  first  of  Eq.  (AlOa)  .  The  modal 
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character  of  7Jj{n)f  for  the  complex  depicted  in  Fig.  3,  as  a 

function  of  {a>/a}J  in  (Ao)/6)J  and  for  a  modal  overlap  parameter 

h  is  exemplified  in  Fig.  9a.  The  modal  character  in 

^nrr  —  ^5  '  ^ 

=  is  clearly  demonstrated  in  this  figure.  It  emerges 
that  were  the  modal  overlap  parameters  (K)  above  unity,  the 

modes  are  no  longer  distinguishable  and  the  values  of  %(«),  for 
discrete  values  of  the  index  (r)  and  a  fixed  value  of  the  index 
(n)  ,  is  reasonably  smooth  as  a  function  of  {col co^  xn  {Aa)l (o^)  .  The 

presence  of  modes  in  as  a  function  of  {col co^)  xn  (Aco! o}„) ,  is 

d 

obscured.  Indeed,  the  non~modal  (smooth)  character  of  TJj{ri),  for 

d 

the  complex  depicted  in  Fig.  3;  namely  for  (1)  and  for  a  modal 
overlap  parameter  b^=b^=2.0,  is  exemplified  in  Fig.  9b.  The 


absence  of  modal 

character 

in  rfj  (1) 

in  Fig. 

9b,  is. 

thus,  clearly 

demonstrated. 

[cf.  Fig. 

9a.  ]  1 

Moreover, 

it  follows  that  for 

modal  overlap 

parameters 

that 

exceed 

unity; 

b^  >  1.0  ,  the 

replacement  of  the  summation  in  Eq.  (AlOa)  by  the  integration  in 
Eq.  (AlOb)  is  naturally  qualified.  What  happens  if  the  modal 

overlap  parameters  are  below  unity;  namely,  when  b^^l.O  ?  It  is 
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argued  that  for  discrete  values  of  the  index  (r)  and  a  fixed 


value  of  the 

index 

(n)  , 

TJj{n,r)  is  not 

a  smooth  summand  and. 

therefore,  in 

this 

case. 

the 

replacement 

of  a  summation  by  an 

integration  is 

not 

naturally 

qualified. 

However,  if  the  index 

(r)  is  made  into  a  continuous  variable,  this  action,  as  the 
second  of  Eq.  (AlOa)  suggests,  is  tantamount  to  planting 

multiplicity  of  modes  in  between  the  originally  discrete  modes. 
These  generated  modes,  as  do  the  initial  modes,  carry  modal 

extents  that  are  largely  determined  by  loss  factors  of  values 

given  by  +  •  Clearly,  the  increase  in  the  modal 

density  due  to  this  transition  from  the  discrete-to-the 
continuous  (r)  renders  T}j{n,r)  smooth  and,  therefore,  qualifying 

the  replacement  of  the  summation  by  integration.  In  this 
procedure  the  artificial  increase  in  the  modal  density  is, 

however,  compensated  by  averaging  over  the  number  of  added  modes, 
as  exemplified  in  Eqs .  (A9a)  -  (Alla).  The  result  of  replacing  a 

summation  by  an  integration  is  commensurate  with  an  artificial 
increase  in  the  modal  overlap  parameter  to  values  that  subdue 

modal  recognitions.  Indeed,  the  a  priori  replacement  of  a 
summation  by  an  integration  results  in  values  of  the  induced  loss 

factor  from  nW)  to  7/W  that  match  these  obtained  by  the  mean- 
value  averaging  proposed  by  Skudrzyk  [14].  [cf.  Eq.  (AlOa).] 
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In  EA:  (Master  in  isolation) 


Fig.  1.  An  externally  force-driven  isolated  master  dynamic 
system. 

rig  =  input  power  generated  by  an  external  force-drive. 

=  stored  energy  ^  N^S^}  6°  =  modal  stored  energy. 

=  number  of  modes  =  A(o)Vg  •  ~  density. 

=  global  mass  =  NJft^  ;  =  modal  mass. 

77^  =  indigenous  loss  factor;  o)  =  frequency  variable. 
A(o)  =  frequency  bandwidth. 


50 


InEA:  (Master) 


(Adjunct) 


Fig.  2.  An  externally  force-driven  master  dynamic  system 
passively  coupled  to  an  adjunct  dynamic  system. 

Subscript  (o)  designates  quantities  and  parameters  that 
pertain  to  the  master  dynamic  system. 

Subscript  (s)  designates  quantities  and  parameters  that 
pertain  to  the  adjunct  dynamic  system. 

Oj  =  the  net  power  that  is  imparted  to  the  adjunct  dynamic 
system  from  the  master  dynamic  system. 

The  vector  k^,  g}  and  its  normalized  form  {jTlc,  (X^,  g"} 

describe  the  coupling  coefficients  {mass,  stiffness  and 
gyroscopic}  between  the  two  dynamic  systems.  [cf.  Fig.  1.] 
rig  =  input  power  generated  by  an  external  force-drive. 

Eg  =  stored  energy  in  the  master  dynamic  system. 

=  stored  energy  in  the  adjunct  dynamic  system. 

=  number  of  modes  in  adjunct  dynamic  system  =  A(<»)V,  ; 

Vy  =  modal  density  in  adjunct  dynamic  system. 

Tj^  =  indigenous  loss  factor  in  the  adjunct  dynamic  system. 

=  global  mass  of  the  adjunct  dynamic  system  =  ; 

=  modal  mass  in  the  adjunct  dynamic  system. 
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Oscillators 
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Fig.  3.  A  Sketch  of  a  specific  complex  dynamic  system 
comprising  a  master  dynamic  system  consisting  of  a  single 
harmonic  oscillator  (a  master  oscillator)  and  an  adjunct  dynamic 
system  consisting  of  several  harmonic  (satellite)  oscillators. 
The  satellite  oscillators  are  identified  by  the  index  (r) .  The 
satellite  oscillators  are  uncoupled  to  each  other.  Also,  only 
the  master  oscillator  is  externally  force-driven. 
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X°(r: 


Fig.  4a.  The  normalized  resonance  frequency  distribution  X°(r)  of 
the  harmonic  oscillators  in  the  adjunct  dynamic  system  (the 
satellite  oscillators)  as  a  function  of  the  index  (r) .  X°(r)  is 


as  stated  in  Eq.  (29a)  and  iV  =27.  0  Discrete, 


Continuous . 


Fig.  4b.  The  localized  loss  factor  //(r)  of  the  harmonic 
oscillator  in  the  adjunct  dynamic  system  (the  satellite 
oscillators)  as  a  function  of  the  index  (r) .  as  stated  in  Eq 

(29b),  N^=rj  and 

O  Discrete,  -  Continuous  for  Z?  =  0.1, 

X  Discrete,  -  Continuous  for  b  =  2.0, 

□  Discrete, -  Continuous  for  b  —  lO. 
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Fig.  5. 


The  induced  loss  factor  Tjj{(o)  as  a  function  of  the 

normalized  frequency  {co/co^  for  N^=27,  {m^I mJ)  =  0.1  and  for  three 
values  of  the  coupling  factor  C{q))  ;  C(£y)  =  1.0  (solid  curve), 

C{co)=3xl0  (dash  curve),  and  C(©)  =  10  (dash-dot  curve).  The 

frequency  bandwidth  [A(6>)/6>J  =  (//2)e0.3. 


Fig.  6.  The  indigenous  loss  factor  7]^^)  in  the  adjunct  dynamic 
system  as  a  function  of  the  normalized  frequency  for 

=27,(M^/Mp)  =  0.1  and  for  three  values  of  the  modal  overlap 
parameter  b^ifo)}  b^{a:^=0.\  (solid),  b^{o)^  =  '1.0  (dash)  and  Z?^(£y)  =  10.0 
(dash-dot).  The  frequency  bandwidth  \_is{ci))l CDq^  =  {y /'2)  =  0.3. 


Fig.  7a.  The  modal  coupling  strengths  as  a  function  of  the 

normalized  frequency  {oloj^)  for  a  coupling  factor 
C(<y)=1.0,  N^=27,  {m^  I  M^)  =  0.\  and  for  three  values  of  the  modal 
overlap  parameter  b^{(a)=0A  (solid  curve),  b^{o))  =  2.0  (dash 

curve)  and  Z?y(6?)  =  10.0  (dash-dot  curve)  . 
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Fig.  7b.  The  modal  coupling  strengths  ^  function  of  the 

normalized  frequency  for  a  coupling  factor 

C(o)  =  3xl0~^,  N^=n,  (MJ Mg)  =  0.\  and  for  three  values  of  the  modal 
overlap  parameter  =  0.1  (solid  curve) ,  Z?^(<y)  =  2.0  (dash 

curve)  and  Z?y(ffl)=10.0  (dash-dot  curve)  . 
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Fig.  7c.  The  modal  coupling  strengths  ®  function  of  the 

normalized  frequency  {o)/co^)  for  a  coupling  factor 

-“3 

C(6>)  =  10  ,  N^=21,  {mJ M^  =  Q.\  and  for  three  values  of  the  modal 
overlap  parameter  Z?^(o)  =  0.1  (solid  curve),  Z?^(6))  =  2.0  (dash 

curve)  and  Z?^(o)  =  10.0  (dash-dot  curve)  . 


Fig.  8.  The  minimum  value  of  the  modal  overlap  parameter 

as  a  function  of  the  normalized  frequency  (cola,,),'  for 

N^=27  and  {m^/M„)  =  0.1  and  for  three  values  of  the  coupling 

factor  C{a);  C(fi>)=l  (dash-dot),  C{a)  =  3x\0  (dash)  and  C(£o)  =  10 

(solid) .  The  statistical  energy  analysis  is  sufficiently  valid 
for  corresponding  values  of  the  modal  overlap  parameter  b^(a)  that 
exceed  the  minimum  value  shown  in  this  figure. 
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Fig.  9.  Discrete  induced  loss  factor  Tjj\\)  as  a  function  of 
(colcOg),  for  a  stiffness  control  coupling  form  with 

=  1 .0 [a  =  0.0.]  ,  g  =  mc=0  [Sprung  masses].  The  modal  distribution 
in  the  adjunct  dynamic  system  is  as  stated  in  Eq.  (29a)  and  b(r) 
in  Eq.  (29b}  is  given  by  O. 

d  . 


with  Z?  =  (0.1) 


b.  With  h  =  (2.0) 
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